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f Abstract. We report a measurement of cosmic shear correlations using an effective area of 6.5 sq. deg. of the 

1~ \ VIRMOS deep imaging survey in progress at the Canada-France-Hawaii Telescope. We measured various shear 

correlation functions, the aperture mass statistic and the top- hat smoothed variance of the shear with a detection 
significance exceeding 12 a for each of them. We present results on angular scales from 3 arc-seconds to half a 
degree. The consistency of different statistical measures is demonstrated and confirms the lensing origin of the 
signal through tests that rely on the scalar nature of the gravitational potential. For Cold Dark Matter models 
, £^ we find os^o' 6 — 0.431q'o5 a t the 95% confidence level. The measurement over almost three decades of scale 

allows to discuss the effect of the shape of the power spectrum on the cosmological parameter estimation. The 
degeneracy on ag — flo can be broken if priors on the shape of the linear power spectrum (that can be parameterized 
by F) are assumed. For instance, with F = 0.21 and at the 95% confidence level, we obtain 0.6 < as < 1-1 and 
0.2 < flo < 0.5 for open models, and a% > 0.65 and Qo < 0.4 for fiat (A-CDM) models. From the tangential/radial 
modes decomposition we can set an upper limit on the intrinsic shape alignment, which was recently suggested 
as a possible contribution to the lensing signal. Within the error bars, there is no detection of intrinsic shape 
alignment for scales larger than 1'. 
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1. Introduction 

Cosmological gravitational lensing produced by 
large-scale structure (or cosmic shear) has been 
advocated as a powerful tool to probe the mass 
distribution in the universe (see the reviews from 



Mellier 1999, Bartclmann & Schneider 2001 and refer 
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year (Van Waerbeke et al. 2000, 


Bacon et al. 2000 


Kaiser et al. 200C 


Wittman et al. 2000 


Maoli et al. 2001 




Rhodes et al. 2001) 


confirmed that 



the amplitude and the shape of the signal are compatible 
with theoretical expectations, although the data sets 
were not large enough to place useful constraints on 
cosmological models. Maoli et al. 2001 combined the 



results from different groups to obtain constraints on the 
power spectrum normalization as and the mean density 
of the universe f2o : Their result is in agreement with the 
cluster abundance constraints, but they were not yet able 
to break the degeneracy between as and £!q. 
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The physical interpretation of the weak lensing sig- 
nal can be made more securely using detections of cos- 
mic shear from different statistics and angular scales 
on the same data set (as in Van Wacrbckc et al. 200Cj[ ) . 
Unfortunately, their joint detection of the variance and 
the correlation function using the same data was not fully 
conclusive: the sample was too small to enable a signif- 
icant detection of the cosmic shear from variances with 
different weighting schemes and 2-points statistics over a 
wide range of scales. The use of independent approaches 
is nevertheless necessary and it is a crucial step to val- 
idate the reliability of cosmic shear, to check the con- 
sistency of the measurements against theoretical predic- 
tions and to understand the residual systematics. A rel- 
evant example is the aperture mass statistic (defined in 
Schneider et al. 1998 ). It is a direct probe of the projected 
mass power spectrum, and it is not sensitive to certain 
type of systematics (like a uniform PSF anisotropy) which 
may corrupt the top-hat smoothed variance, or the shear 
correlation function. Even the shear correlation function 
can be measured in several ways, by splitting the tangen- 
tial and radial modes for instance. 

In this paper we report the measurement of the top-hat 
smoothed variance, the aperture mass, the shear correla- 
tion function, and the tangential and radial shear corre- 
lation functions on a new homogeneous data set covering 
an effective area of 6.5 square-degrees (deg 2 ). The depth 
and the field of view are well suited for a comprehensive 
analysis using various statistics. We show that the ampli- 
tude of residual systematics is very low compared to the 
signal and discuss the consistency of these measurements 
against the predictions of cosmological models. 

We also discuss alternative interpretations. It 
has been suggested recently that intrinsic align- 
ments of galaxies caused by ti dal fields could con- 
tr ibute to the lensing signal (Croft fc Mctzlcr 2000 



Heavens, Refregier fc Heymans 2000[ |Catclan et al. 200C 



Crittenden et al. 2000a| , |Crittendcn et al. 2000b| ). This 



type of systematic is problematic because its signature 
on different 2-points statistics mimics the lensing effect. 
A mode decomposition in electric and magnetic types 
(or E and B modes), similar to what is performed 
for the polarization analysis in the Cosmic Microwave 
Background, can separate lensing from intrinsic alignment 
(see |Crittenden et al. 2000al prittenden et al. 2000b| ). 



The E and B mode analysis is the subject of a forth- 
coming paper; the aperture mass statistic presented in 
this paper is a similar analysis to the E and B mode 
decomposition, and allows us to put an upper limit on the 
contamination of our survey by the intrinsic alignments. 

This paper is organized as follow: Section 2 de- 
scribes our data set, and highlights the differences 
in th e data preprocessing f rom our previous analy- 



sis ( |Van Waerbeke et al. 2000| ). The measurement of the 
shear from this imaging data is discussed in Section 3. 
Section 4 summarizes the theoretical aspects of the differ- 
ent quantities we measure, and lists the statistical estima- 
tors used. The results and comparison to a few standard 



cosmological models are shown in Section 5. In Section 6 
we perform a maximum likelihood analysis of cosmologi- 
cal models in the (Oo, as) parameter space. The results on 
very small scales are shown separately in Section 7, and 
we conclude in Section 8. 



2. The data set 

The DESCART weak lensing project]] is a theoretical and 
observational program for cosmological weak lensing in- 
vestigations. The cosmic shear survey carried out by the 
DESCART team uses the CFH12K data jointly with the 
VIRMOS survey^ to produce a large homogeneous photo- 
metric sample which will eventually contain a catalog of 
galaxies with redshifts a s well as the projecte d mass den- 
sit y over the whole field ( Le Fevre et al 2001 ) . In contrast 
to 



Van Waerbeke et al. 200C , the new sample presented in 



this work only uses I-band data taken with the CFH12K 
camera and is therefore more homogeneous. It is worth 
noting that our new CFH12K sample only uses half of the 
data of the previous one. A comparison of the results will 
also permit to check the consistency and the robustness 
of the cosmic shear analysis. 

The CFH12K data was obtained during dark nights 
in May 1999, November 1999 and April 2000 follow- 
ing the standard observation procedure described in Van 
Waerbeke et al (2000). The fields are spread over 4 in- 
dependent 2x2 deg 2 areas of the sky identified as F02, 
F10, F14 and F22. Each field is a compact mosaic of 16 
CFH12K pointings named P[n] n= i_i6. Once the survey is 
completed each of them will cover 4 deg 2 . Currently, of 
the final 16 deg 2 , only 8.38 deg 2 is available for the anal- 
ysis - most of the pointings are located in three different 
fields (F02, F10, F14 listed in Table 0). This total field of 
view gets significantly reduced by the masking and selec- 
tion procedures described below. A summary of the data 
set characteristics are listed in Table |l|. 

The data reduction was done at the TERAPIX data 
centerf] More than 1.5 Tbytes of data were processed in 
order to produce the final s tacked images. The reduct ion 
procedure is the same as in Van Waerbeke et al. 200C| , so 
we refer to this paper for the details. However, in order 
to improve the image quality prior to correction for the 
PSF anisotropy and to get a better signal-to-noise ratio 
on a larger angular scale than in our previous work, all 
CFH12K images were co-added after astrometric correc- 
tions. 

The astrometric calibration and the co-addition were 
done using the MSCRED package in IRAF. Some tasks 
have been modified in order to allow a fully automatic 
usage of the package. For each pointing, we first started 
with the images in the I band. An astrometric solution was 
first found for one set of exposures in the dither sequence 
using the USNO-A 2.0 as reference, which provides the 



http:/ /terapix. iap.fr/Descart 
http:/ /www. astrsp-mrs.fr/virmos/ 



http:/ /terapix. iap.fr 
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Table 1. List of the fields. All observations were done in I band with the CFH12K camera (Cuillandre et al. 200C). The 
number following the F denotes the field name, and the number following the P denotes the pointing name within the 
field. The geometry of the survey is detailed in http://terapix.iap.fr/Dcscart/. The image quality has been measured 
on each stacked image from a standard fitting of a Moffat profile. 



Target 


Used area 


Exp. 


time 


Period 


Image quality 


F02P1 


980 arcmin 2 


9390 


sec. 


Nov. 1999 


0.75" 


F02P2 


1078 arcmin 2 


7200 


sec. 


Nov. 1999 


0.90" 


F02P3 


980 arcmin 2 


7200 


sec. 


Nov. 1999 


0.90" 


F02P4 


1078 arcmin 2 


7200 


sec. 


Nov. 1999 


0.80" 


F10P1 


882 arcmin 2 


3600 


sec. 


May 1999 


0.65" 


F10P2 


882 arcmin 2 


3600 


sec. 


May 1999 


0.75" 


F10P3 


490 arcmin 2 


3600 


sec. 


May 1999 


0.75" 


F10P4 


882 arcmin 2 


3600 


sec. 


May 1999 


0.65" 


F10P5 


882 arcmin 2 


3600 


sec. 


May 1999 


0.75" 


F10P7 


1176 arcmin 2 


3600 


sec. 


Apr. 2000 


0.75" 


F10P8 


1176 arcmin 2 


3600 


sec. 


Apr. 2000 


0.70" 


F10P9 


98 arcmin 2 


3600 


sec. 


Apr. 2000 


0.65" 


F10P10 


784 arcmin 2 


3600 


sec. 


Nov. 1999 


0.80" 


F10P11 


294 arcmin 2 


3600 


sec. 


Nov. 1999/Apr. 2000 


0.90" 


F10P12 


1176 arcmin 2 


3600 


sec. 


Apr. 2000 


0.80" 


F10P15 


686 arcmin 2 


3600 


sec. 


Apr. 2000 


0.85" 


F14P1 


882 arcmin 2 


3600 


sec. 


May 1999 


0.80" 


F14P2 


882 arcmin 2 


3600 


sec. 


May 1999 


0.85" 


F14P3 


686 arcmin 2 


3600 


sec. 


May 1999 


0.75" 


F14P4 


1078 arcmin 2 


3600 


sec. 


May 1999 


0.75" 


F14P5 


980 arcmin 2 


3600 


sec. 


May 1999 


0.70" 


F14P6 


686 arcmin 2 


3600 


sec. 


May 1999 


0.80" 


F14P7 


686 arcmin 2 


3600 


sec. 


May 1999 


0.70" 


F14P8 


882 arcmin 2 


3600 


sec. 


May 1999 


0.85" 


F14P9 


1078 arcmin 2 


3600 


sec. 


Apr. 2000 


0.75" 


F14P10 


784 arcmin 2 


3600 


sec. 


May 1999 


0.85" 


F14P11 


882 arcmin 2 


3600 


sec. 


Apr. 2000 


0.80" 


F14P12 


784 arcmin 2 


3600 


sec. 


Apr. 2000 


0.80" 


F14P13 


882 arcmin 2 


3600 


sec. 


Apr. 2000 


0.85" 


F14P14 


882 arcmin 2 


3600 


sec. 


May 1999 


1.0" 


F14P15 


882 arcmin 2 


3600 


sec. 


Apr. 2000 


0.90" 


F14P16 


1176 arcmin 2 


2880 


sec. 


Apr. 2000 


0.65" 


F22P3 


686 arcmin 2 


3600 


sec. 


May 1999 


0.75" 


F22P4 


980 arcmin 2 


3600 


sec. 


Nov. 1999 


0.75" 


F22P6 


588 arcmin 2 


3600 


sec. 


Apr. 2000 


0.80" 


F22P11 


294 arcmin 2 


2880 


sec. 


Apr. 2000 


0.75" 



position of ~ xlO 8 sources with an RMS accuracy of 0.3 
arcsec. The astrometric solution was then transferred to 
the other exposures in the sequence. All object catalogs 



were obtained using SExtractor ( Bcrtin fc Arnouts 199C| )| | 
and a linear correction to the world coordinate system 
was computed with respect to the initial set. Finally, all 
images were resampled using a bi-cubic interpolation and 
then stacked together. 

At this stage, each stacked image was inspected 
by eye and all areas which may potentially influ- 
ence the later lensing an alysis signal were mask ed (see 
Van Waerbeke et al. 2000| and [Maori et al. 200l[ ). Since 
we adopted conservative masks, this process had a dra- 
matic impact on the field of view: we lost 20% of the total 
area and ended up with a usable area of 6.5 deg 2 . 



The photometric calibrations were done using stan- 
dard stars from the Landolt's catalog ( |Landolt 1992 ) cov- 
ering a broad sample of magnitude and colors. A full de- 
scription of the photometric procedure is beyond the scope 
of this work and will be discussed elsewhere (Le Fevre et 
al, in preparation). In summary, we used the SA110 and 
SA101 star fields to measure the zero-points and color 
equations of each run. From these calibrations, we pro- 
duced the magnitude histograms of each field in order 
to find out the cut off and a rough limiting magnitude. 
Although few fields have exposure time significantly larger 
than 1 hour, the depth of the sample is reasonably stable 
from field to field and reaches Iab = 24.5. Up to this 
magnitude, 1.2 million galaxies were detected over the to- 
tal area of 8.4 deg 2 . 



4 ftp:/ /ftp. iap.fr/pub/from_users/bertin/sextractor/ 
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3. Shear measurements 

3.1. Shape measurement 

The details of our shape measurement procedure 
and Point Spread Function (hereafter PSF) correction 
have been extensively described in two previous pa- 
pers (|Van Waerbeke et al. 2000[ [Maoli et al. 200l|), and 



tested against numerical simulations (Erbcn et al. 2001). 
Therefore we will not reproduce these details here, 
but only give a short overview of the procedure. The 
s hape measuremen t pipeline uses the IMCAT software 
( Kaiser et al. 1995| )P1 combined with the SExtractor pack- 
age. The different steps in the procedure are as follows: 



Object detection with Sextractor. 



arc 



The shape parameters defined in Kaiser et al. 1995 
calculated using IMCAT. 

Stars are identified in the stellar branch of the size- 
magnitude diagram. Stars brighter than 1 magni- 
tude below the saturation level are excluded. Objects 
smaller than the PSF size are discarded as galaxy can- 
didates (because a shape correction below the PSF size 
is meaningless). 

The PSF is measured from the stars, and interpolated 
continuously over the CCD's using a third order poly- 
nomial. 

Galaxy shapes are corrected using the scheme devel- 
oped in Kaiser et al. 1995 m odified and adapt ed to 
our problem as described in Erben et al. 2001, The 
shape correction is a two-step process: first we re- 
move the anisotropic contribution of the PSF, then 
the isotropic contribution is suppressed according to 
Luppino & Kaiser 1997 



— A weight w is calculated for each galaxy, which de- 
pends on the level of noise in the shape correction (see 
Eq.(7) in [Van Waerbeke et al. 2000| ). 

— For each galaxy pair with members closer than 15 pix- 
els (3 arcsec), one member is removed, in order to 
avoid the problem of overlapping isophotes reported 
in |Van Waerbeke et al. 20"0C| . 

— Each CCD is visualized by eye, and the bad areas are 
masked (star spikes and ghost images, blank lines or 
columns, fringe residuals). After the whole process of 
cleaning and object selection, 420,000 galaxies were 
effectively used for the weak lensing analysis. 

The raw ellipticity e of a galaxy is measured from the 
second moments Iij of the surface brightness f(9): 



e = 



hi - I: 



2L 



Tr( I ! ' 7>(ij) ' 



in = / d A 9W(e)9 i e j f(e).{i) 



The window function W{6) suppresses the noise at large 
distances from the object center. The procedure described 
above gives a corrected galaxy ellipticity e gal calculated 



http: / /www. ifa.hawaii.edu/~kaiser/ 
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Fig. 1. Top and third from the top panels: averaged ef al 
and ef al component versus the ellipticity of the PSF at the 
galaxy location, without the anisotropic PSF correction. 



Second from top and bottom panels: averaged ef al and 
e| al including the anisotropic PSF correction. 



from the e's. According to Kaiser et al. 1995, the ensem- 



ble average of e gal is equal to the shear 7 at the galaxy lo- 
cation. Figure |l] shows the level of systematics in e gal with 
and without the anisotropic PSF correction. After the PSF 
correction, the average galaxy ellipticity is bounded be- 
tween —0.005 and 0.005, and the variance of the resid- 
ual systematics is less than ~ 10 -5 . As we shall see later 
this is much less than the measured signal. As quoted in 
Van Waerbeke et al. 2000 , the galaxy ellipticities show a 



small offset (—0.008,-0.003), which has been corrected 
in this figure (the origin of this offset is still unclear). 
However it is worth to mention that the aperture mass 
is not sensitive to this offset. 



4. Statistical measures of shear correlations 

4.1. Theory 

We summarize the different statistics we shall measure, 
and how they depend on cosmological models. We con- 
centrate on 2-points statistics and variances, since higher 
order moments are more difficult to measure, and will be 
addressed in a forthcoming paper. 

Let us assume a source redshift distribution parame- 
terized as: 



kindly made available b y Nick Kaiser at n[Z s ) 



(3 



z r 



— exp 



Zs\ 

zo) 



(2) 
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where the parameters (zo,a,/3) = (0.8,2,1.5), which is 
consistent with a limiting magnitude Iab — 24.5 given 
by Cohen et al. 200C (it corresponds to a mean redshift of 
1.2). We define the power spectrum of the convergence as 
(following the notation in Schneider et al. 1998): 



P K (k) 



4 : U 



P 3 D 



k 



Ik{w) 



, w 



dw n(w ) — , 
Tk{W) 



(3) 



where /if(u>) is the comoving angular diameter distance 
out to a distance w (wu is the horizon distance), and 
n(w(z)) is the redshift distribution of the sources given in 
Eq.(g). fVo(fc) is the non- linear mass power spectrum, and 
k is the 2-dimensional wave vector perpendicular to the 
line-of-sight. For a top- hat smoothing window of radius 
9 C , the variance is: 



(7 2 



2 

^2 



dk 



P K {k)[J x {k6 c )] 2 , 



(4) 



where J\ is the first Bessel function of the first kind. 
The aperture mass M ap was introduced in 



Kaiser et al. 1994 



Map = 



d 2 e K {e) u(0), 



(5) 



where k(8) is the convergence field, and U(9) is a com- 
pensated filter (i.e. with zero mean). [Schneider et al. 1998 
applied this statistic to the cosmic shear measurements. 
They showed that the aperture mass variance is related to 
the convergence power spectrum by: 



(Ml p ) 



288 



dk 
F 



P K (k)[J 4 {k9 c )f 



(6) 



(M ap ) can be calculated directly from the shear 7 without 
the need for a mass reconstruction. 

For each galaxy, we define the tangential and radial shear 
components (74 and 7 r ) with respect to the center of the 
aperture: 



7t = -7i cos(20) - 72 sin(2(/>) 
lr = -72cos(20) +7isin(20), 



(7) 



where <f> is the position angle between the x-axis and the 
line connecting the aperture center to the galaxy. It is 
then easy to show that the aperture mass is related to the 
tangential shear by: 



Map = 



d 2 7t (0) Q(9), 



(8) 



<0c 



where the filter Q(9) is given from U(6): 



U(d) 



If 7t is replaced by 7 r in Eq.(|J), then the lensing sig- 
nal vanishes, due to the curl- free property of the shear 
field (Kaiser et al. 1994) ^. This remarkable property con- 
stitutes a test of the lensing origin of the signal. The 
change from 7t to j r can simply be accomplished just by 
rotating the galaxies by 45 degrees in the aperture (i.e. 
changing a curl- free field to a pure curl field). Hereafter 
we call the M ap statistic measured with the 45 degrees 
rotated galaxies the i?-mode (R for radial mode), and 
(M'j_) the corresponding variance. It is interesting to note 
that the i?-mode is not expected to vanish if the mea- 
s ured signal is due to intrinsic alignments of galaxies 
( ICrittenden et al. 2000b|) . Therefore it can be used to con- 
strain the amount of residual systematics as well as the 
degree of the intrinsic alignment of galaxies. 
From the shear 7 and its projections defined in Eq.(0) 
we can also define various galaxy pairwise correlation 
functions related to the convergence power spectrum. 
Note that the tangential and radial shear projections in 
what follows are performed using the relative location 
vector of the pair members, not from an aperture cen- 
ter. The following correlation functions can be defined 
( [Miralda-Escude 199l| , [Kaiser 1992j ): 



1 f°° 
{ll)9 = ^J o dk P - 



(k)J (k9), 



(ltlt)t 



2vr 



<7r7r>* = ^ 



dk P K (k)[J (k9) + Ji(k8)], 



dk P K (k)[J (k9) - J 4 (k9)}, 



(10) 



(11) 



(12) 



where 9 is the pair separation angle. The cross-correlation 
(7t7r)e is expected to vanish for parity reasons (there is 
no preferred orientation on average). 

It is easy to see that the Eqs.(|j|,|l|,|ll|,|l2|) are differ- 
ent ways to measure the same quantity, that is the con- 
vergence power spectrum P K (k). Ultimately the goal is to 
deproject P K (k) in order to reconstruct the 3D mass power 
spectrum from Eq.(||), but this is beyond the scope of this 
paper. Here we restrict our analysis to a joint detection 
of these statistics, and show that they are consistent with 
the gravitational lensing hypothesis. We will also examine 
the constraints on the power spectrum normalization as 
and the mean density of the universe Qq ■ 

4.2. Estimators 

Let us now define the estimators we used to measure the 



quantities given in Eqs.flfflJIOWlM. 

The variance of the shear is simply obtained by a cell 
averaging of the squared shear j 2 (0i) over the cell index 



6 Curl modes are p roduced by non-linear le nsing effects, but 
v / these are very small (Bernardeau et al. 1997). 
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An unbiased estimate of the squared shear for the cell i where ef al is the tangential galaxy ellipticity, and Q is 



is: 



2 iVi 



et{0k)et{0i) 

Wkwi 

k^l 



(13) 



where Wk is the weight for the galaxy fc, and Ni is the 
number of galaxies in the cell i. The cell averaging over 
the survey is then an unbiased estimate of the shear vari- 
ance (7 2 ). However, due to the presence of masked areas 
(mentioned in Section 4.1), some cells may have a very 
low number of galaxies compared to others. Instead of ap- 
plying an arbitrary sharp cut off on the fraction of the 
apertures filled with masks (as it was in previous works) 
we decided to keep all the cells, and to weight each of them 
with the squared sum of the galaxy weights located in the 
cell. The cell averaging is now defined as: 



E[l 2 



E 

cells 



\fe=i / 



E 

cells 




(14) 



where i identifies the cell. One potential problem with 
this procedure is that the sum of the weights is related to 
the number of objects in the aperture, which is affected by 
magnification bias, and therefore correlated with the shear 
signal measured in the same aperture. Fortunately the 
first non-vanishing contribution of this weighting scheme 
is a third order effect (of order 1%), and is therefore 
negligible^. The advantage is that we can use all cells 
without wondering about their filling factor, and it natu- 
raly down-weights the cells which contain a large fraction 
of poorly determined galaxy ellipticities. The weighting 
scheme Eq.([l4]) has been tested against numerical simu- 
lation, using a simulated survey with exactly the same 
survey geometry as our data: it gave unbiased measures 
of the lensing signal applied to the galaxies. 

The M ap statistic is calculated from a similar estima- 
tor, although the smoothing window is no longer a top-hat 
but the Q function defined in Eq. (^J) . An unbiased estimate 
of in the cell i is: 



E[Mt v {0i)] = 



Ni 

E 

k^l 



w k wi e? al (0 fc )ef al (0 ; ) Q(# fc )Q(0,) 



N 

E 



(15) 



w k wi 



7 Moreover the slope of number counts in our I-band is 
0.3, which makes the magnification effect very small (see 



Moessner, et al. 1998 for an application of the effect to the an- 
gular correlation function). 



given by (see [Bclmcidcr ct al. 1998] ): 



6/0 



Q(0) = - hr 

7T V r 



1- ~ 



(16) 



The estimation of (M^ p ) over the survey is then given 
by the same expression as in Eq.(^4|), with E[j 2 (9i)} re- 
placed by E[M^ p (6i)]. We emphasize that the this filter 
probes effective scale s c /5, and not 9 C (see Figure 2 in 
[Schneider et al. 1998 ). Therefore we have to be careful 
when comparing the signal at different scales between dif- 
ferent estimators. 

The shear correlation function (77) g at separation is 
obtained by identifying all the pairs of galaxies falling in 
the separation interval [8 — dd, 8 + dd] , and calculating the 
pairwise shear correlation: 



2 £ ™ ^m&ipi 
e[ 1T a = Y.— 



pairs 



(17) 



The tangential and radial correlation functions {jt"ft)e and 
(lrlr)e are measured also from Eq.(|l7j) by replacing e gal 
with ef al and ef al respectively and dropping the sum over 
a. It is worth noting that the estimators given here are 
independent of the angular correlation properties of the 
source galaxies. 

5. Results and comparison to cosmological models 

In this section we present our measurements of the 2-point 
correlations of the shear using the different estimators de- 
fined above. Figures || to show the results for the differ- 
ent estimators: the variance in Figure 0, the mass aper- 
ture statistic in Figure |j| the shear correlation function in 
Figure 4, the radial and tangential shear correlations in 
Figure 5, and the cross-correlation of the radial and tan- 
gential shear in Figure ^. Along with the measurements we 
show the predictions of three cosmological models which 
are representative of an open model, a flat model with cos- 
mological constant, and an Einstein-de Sitter model. The 
amplitude of mass fluctuations in these models is normal- 
ized to the abundance of galaxy clusters. The three models 
are characterized by the values of Qo, A and as as follows: 

— short-dashed line: Sl = 0.3, A = 0, erg — 0.9 

— solid line: fl = 0.3, A = 0.7, a 8 = 0.9 

— long-dashed line: flo = 1, A = Q, erg = 0.6 

The power spectrum is taken to be a cold dark matter 
(CDM) power spectrum with shape parameter T — 0.21. 
The predictions for shear correlations are computed us- 
ing the non-linear evolution o f the power spectrum us- 
ing the Peacock fc Dodds 1996 fitting formula. The source 
redshift distribution follows Eq.(||) with (z ,a,f3) = 
(0.8, 2, 1.5), which corresponds to a mean redshift of 1.2. 
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4 5 6 7 8 

0[arcmin] 

Fig. 2. Top-hat smoothed variance of the shear (points 
with error bars). The three models correspond to 
(f7 ,A,cr 8 ) = (0.3, 0,0.9), (0.3, 0.7, 0.9), (1,0, 0.6) for the 
short-dashed, solid and long-dashed lines respectively. The 
power spectrum is a CDM-model with T = 0.21. The error 
bars correspond to the dispersion of the variance measured 
from 200 realizations of the data set with randomized ori- 
entations of the galaxy ellipticities. 



It is reassuring that the different statistics agree with 
each other in their comparison with the model predic- 
tions. These statistics weight the data in different ways 
and are susceptible to different kinds of systematic errors. 
The consistency of all the 2-point estimators suggests that 
the level of systematics in the data is low compared to the 
signal. A further test for systematics is provided by mea- 
suring the cross-correlation function (jtl^e, which should 
be zero for a signal due to gravitational lensing. Figure 
shows the measured cross-correlation function, which is 
indeed consistent with zero on all scales. The figure also 
shows the cross-correlation obtained when the anisotropic 
contamination of the PSF is not corrected - clearly such 
a correction is crucial in measuring the lensing signal. 

The lower panel of Figure || shows the R-mode of the 
mass aperture statistic. As this statistic uses a compen- 
sated filter, the scale beyond which the measured R-mode 
is consistent with zero (5' on the plot) corresponds to an 
effective angular scale 9 ~ 1'. This places an upper limit on 
measured shear correlations due to the intrinsic alignment 
of galaxies, given the redshift distribution of the sources. 
The vanishing of (M±) for effective angular scales larger 
than 1' strongly supports our conclusion that the level 
of residual systematics is low: this is a very hard test to 
pass, as it means that the signal is produced by a pure 
scalar field, which need not be the case for systematics. 




1 ! 



5 6 7 8 

\ arc ; " 



10 11 12 13 



Fig. 3. The aperture mass statistic for the same models 
as in Figure ||. The lower panel plots the R-mode, obtained 
by making a 45 degree rotation as described in the text. 
There is no significant detection for 9 > 5 arcmin (corre- 
sponding to an effective angular scale of 1', as discussed 
in the text), which shows the low level of contamination 
by galaxy intrinsic alignment and/or residual systematics. 



We checked that M\ is Gaussian distributed with a zero 
average all over the survey, as what we would expect from 
a pure noise realisation. For scales below 5' on the plot, 
the i?-mode is not consistent with zero at the 2-a level. 
Since the cross-correlation (j t j r ) e is consistent with zero 
at this scale, the source of the i?-mode is probably not a 
residual system atic. It might be due to t he effect of intrin- 
sic alignments ( Crittenden et al. 2000b| ), but it is difficult 
to be sure without further tests. 

The error bars shown in Figures || to [?] are calculated 
from a measurement of the different statistics in 200 re- 
alizations of the data set, with randomized orientations 
of the galaxies. We measured the sample variance from 
ray-tracing simulations (Jain et al. 2000) and find that it 
is smaller than 20% of the noise error bars shown here 



(see Van Waerbeke et al. 1999 where the sample variance 
has been calculated for surveys with similar geometry), 
therefore we have not included it in our figures. Figure |^ 
shows an estimate of the sample variance for the r.m.s. 
shear using a compact 6.5 sq. degree ray-tracing simula- 
tion (Jain et al. 2000). This figure shows that the sample 
variance is about order of magnitude smaller than the sig- 
nal for the range of scales of interest. Hence our errors 
are not dominated by sample variance, as was the case in 
the first detections of cosmic shear. As the probed angular 
scales approach the size of the fields (which is ~ 30' with 
the CFH12K camera) the sample variance becomes larger. 
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Fig. 4. Shear correlation function (77)0. The models are 
the same as in Figure [|. The lower panel uses a log-scale 
for the x-axis to highlight the small scale details. 
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Fig. 6. Shear r.m.s. (7 2 ) 1 / 2 (solid line) measured in a ray- 



tracing simulation (Jain ct al. 2000) for the open On = 0.3 
model. The dashed line is the sample variance of the shear 
r.m.s. measured from 7 different realisations of the mass 
distribution for a survey of 6.5 sq. degrees. 
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Fig. 5. Top panel: tangential shear correlation function 
{itltje- Bottom panel: radial shear correlation function 
(jrjr)e- The models are the same as in Figure ^. 



This could be responsible for the small fluctuations in the 
measured correlations in Figures |J and || for scales larger 
than 24'. 
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Fig. 7. Shear cross-correlation function {'yt'Jr)e- The sig- 
nal should vanish if the data are not contaminated by 
systematics. As a comparison, the open circles show the 
same cross-correlation function computed from the galaxy 
ellipticities where the anisotropic correction of the PSF 
has been skipped. 

6. Cosmological constraints 



As noted elsew here (e.g. Bernardeau et al. 1997 



Jain fc Seljak 1997| ), the parameters that dominate 
the 2-point shear statistics are the power spectrum nor- 
malization as and the mean density f2o- We investigate 
below how the statistics measured in Figures [| to [| 
constrain these parameters. Our parameter estimates 
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below rely on some simplifying assumptions; a more 
detailed analysis over a wider space of parameters will be 
presented elsewhere. 

We assume that the data follow Gaussian statistics and 
neglect sample variance since it is a very small contribu- 
tor to the noise for our survey, as discussed above. We 
compute the likelihood function C: 



1 



(2tt)™/ 2 |S| 



1/2 



exp 



1 



(d-s) T S- 1 (d-s) 



(18) 



where d and s are the data and model vectors respectively, 
and S := ((d — s) T (d — s)) is the noise correlation ma- 
trix. S was computed for the different statistics from 200 
random realizations of the survey, therefore effects associ- 
ated with the survey geometry are included in our noise 
matrix. The model s was computed for a grid of cosmologi- 
cal models which covers Hq s [0, 1] and erg <E [0.2, 1.8] with 
a zero cosmological constant. The prior is chosen to be flat 
over this grid, and zero outside. We also fixed T — 0.21 
and used the redshift distribution of Eq.(||). We discuss 
below the impact of this choice of priors. 

Figure || (bottom panel) shows that for effective scales 
smaller than 1' there is a non-vanishing R-mode which 
could come either from a residual systematic, or from 
an intrinsic alignment effect. Therefore it is safer to ex- 
clude this part from the likelihood calculation. Thus for 
the top- hat variance, we excluded the point at 1', for the 
correlation functions the points below 2', and for the M ap 
statistic the points below 5'. For the correlation function, 
we also excluded the points at scales larger than 20' be- 
cause of the small fluctuations in the measured correla- 
tions. The constraints on the cosmological parameters are 
not significantly affected whether these large scale points 
are excluded or not. 

Figures || to |lj show the (flo, erg) constraints for each 
of the statistics shown in Figures || to ||. The contours 
show the 99.9%, 95.0% and 68.0% confidence levels. The 
agreement between the contours is excellent, though the 
M ap statistic and the radial correlation function do not 
give as tight constraints as the other statistics. The corre- 
lation function measurements below 2' may be considered 
by using error bars that include a possible systematic bias: 
this is equivalent to adding a systematic covariance ma- 
trix S sys to the noise covariance S matrix in Eq. (|l8|) . The 
new contours computed with the enlarged error bars^j are 
shown in Figure The maximum of the likelihood in 
the variance and correlation function likelihood plots is at 
as — .9 and Sip ~ 0.3. Note that compared to a similar 
plot in Maoli et al. 2001 (Figure 8), here the contours are 
narrower, and are obtained from a homogeneous data set. 
Moreover, the degeneracy between f2 and erg is broken. 

The partial breaking of degeneracy between J7 and 
as was expected from the fully non-linear calculation of 



shear correlations (Jain & Scljak 1997). In the non-linear 
regime the dependence of the 2-points statistics on fio 

8 The enlarged error bars were computed from the estimation 
of our B-mode analysis which will be presented elsewhere 



and as becomes sensitive to angular scale. For example, 
as shown in Jain & Seljak 1997, the shear r.m.s. measures 
as Sip 5 on sca le between 2' — 5', and as fig' 8 on scales 
^10'. Therefore a low fio universe should see a net de- 
crease of shear power at large scale compared to a fio = 1 
universe (for a given shape of the power spectrum) , as ev- 
ident in Figure |[ Note that the aperture mass M ap is still 
degenerate with fi and as (Figure ||) because it probes 
effective scales up to ~ 2.6' only, which is not enough to 
break the degeneracy. 

It seems that the aperture mass (Figure gives a 
slightly larger er 8 for a large f2o compared to the other 
statistics, while they all agree for Qq < 0.7. This could be 
an indication for a low Qq Universe, however in practice, 
the probability contours for the different statistics can- 
not be combined in a straightforward way because they 
are largely redundant. The best strategy here is to con- 
centrate on one particular statistic. We expect the best 
constraints from the shear correlation function (since it 
contains all the information by definition), and therefore 
base our parameter estimates on the likelihood contours 
obtained from it. The contours in the erg — 51o plane in 
Figure [l3] closely follow the curve er 8 cx Qq 6 . This allows 
us to obtain the following measurement of erg fig' 6 (from 
this figure alone): 



a g n° 6 = 0.43 



-0.04(0.06) 
-0.05(0.07)' 



(19) 



where the uncertainties correspond to the 95% (99.9%) 
confidence levels. The result in equation ( |l9| ) is fairly ro- 
bust against different values of T. 

If we fix r = 0.21, we can constrain the two pa- 
rameters separately; we get, at the 95% confidence level: 
0.2 < Sl < 0-5 and 0.6 < erg < 1.1 for open models 
and erg > 0.65 and O < 0.4 for flat (A-CDM) models. 
This result however is sensitive to the prior choosen for 
r. In particular, if we use the relation T = Cloh for a 
cold dark matter model, then some extreme combinations 
of as, Slo and T cannot be ruled out from lensing alone. 
The degeneracy between Oo and erg is broken only if we 
take r to lie in a reasonable interval. Such interval can 
be motivated by galaxy surveys for instance, which give 
0.19 < T < 0.37 at 68% confidence level for the APM 



(Eisentsein & Zaldarriaga 2001). Therefore the separate 
constraints on Oq and as given above require some prior 
assumptions and must be taken with precaution, while 
the constraint on erg S7q 6 is much more robust. The red- 
shift distribution of the sources is likely to be the main 
source of uncertainty in our estimate of equation (|l9|); 
a rough guide is given by the scaling erg Qq 6 cx Zq' 5 



(Jain & Seljak 1997). A more detailed analysis of param- 
eter estimation is left for a later study. 



7. Small scale signal 

Our correlation function measurements extend to much 
smaller scales than shown in the figures above. The limit 
is set only by the fact that we reject one member of all 
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Fig. 8. Likelihood contours in the f2o — 08 plane from 
the top- hat smoothed variance (7 2 ) shown in Figure g. 
The first point in Figure || was not included in the like- 
lihood calculation to avoid the small scale systematic 
shown in Figure || (bottom panel) . The cosmological mod- 
els have A = 0, with a CDM-type power spectrum and 
T = 0.21. The rcdshift of the sources is given by Eq.(Q). 
with (z ,a,(3) = (0.8,2,1.5). The confidence levels are 
(0.68,0.95,0.999). 
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Fig. 9. As in Figure |8J, but using the M ap statistic of 
(Figure ||, top panel) instead of the top-hat variance. The 
first five points in Figure |^ were not included in the likeli- 
hood calculation in order to avoid the small scale system- 
atic shown in Figure || (bottom panel). 



pairs closer than 3 arcsec. Figures |TJ and [l5] show the tan- 
gential, radial and total shear correlation functions. The 
pair separation bins are much smaller than in Figures |] 
and |^, which explains why the error bars are larger. Even 
at the smallest scales, the shear correlation function (77)0 
is consistent with the model predictions. 

The surprising result for the small scale correlations is 
the behavior of the tangential and radial shear correlation 
functions: at scales smaller than 5" we find an increased 
amplitude for (7*74)9, and a negative {^ r ^ r ) g . Though sur- 
prinsing, a negative (•y r 'y r }g is not unphysical: for instance 



b 1.0 




0.5 - 



Fig. 10. Likelihood contours as in Figure 8, but using 
the shear correlation function (77)0 (Figure 1) instead of 
the top-hat variance. The first two points and scales larger 
than 20' in Figure || were not included in the likelihood 
calculation to avoid the contribution from the small scale 
systematic shown in Figure |^ (bottom panel) . 




Fig. 11. As in Figure |3|, but using the tangential shear 
correlation function (7*74)9 (Figure ||) instead of the top- 
hat variance. The first two points and scales larger than 
20' in Figure || were not included in the likelihood calcu- 
lation in order to avoid the contribution from the small 
scale systematic shown in Figure ^ (bottom panel). 



in Kaiser 1992| (Table 1) a shallow mass power spectrum 
(n > —1) implies such an effect. In terms of halo mass 
profile, it corresponds to a projected profile steeper than 
— 1.5. However, regardless of the nature of this signal, it 
is important to note that this is a very small scale ef- 
fect which has no effect on the statistics discussed in pre- 
ceding sections. The contribution of the increased signal 
from (7474)9 to the variance at 1' is less than 1%; more- 
over since (77)9 is not affected at all, the variance is also 
unaffected. As an explicit test, we checked that by remov- 
ing one member of the pairs closer than 7" the measured 
signal in Figures |Hil,il i s unchanged. In a similar cos- 
mic shear analysis using the Red-sequence Cluster Survey 
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Fig. 12. As in Figure ||, but using the radial shear cor- 
relation function (j r ^ r )g (results in Figure ||) instead of 
the top-hat variance. The first two points and scales larger 
than 20' in Figure || were not included in the likelihood 
calculation in order to avoid the contribution from the 
small scale systematic shown in Figure ^ (bottom panel) . 
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Fig. 13. Likelihood contours as in Figure 10, but all the 
points in Figure H on scales smaller than 20' were used. 
In order to account for the small scale systematic shown 
in Figure || (bottom panel) the error bars on the first two 
points were increased to include the systematic amplitude. 



( pladders fc Yee 20pq ) another group finds a similar small 
scale behavior, though at lower statistical significance (H. 
Hoekstra, private communication). 

The cross-correlation (^tlr)e vanishes down to 3", 
therefore no obvious systematic is responsible for this ef- 
fect. The effect is unlikely to be caused by overlapping 
isophotes, or close neighbors effects because {jt)g « 
{ltlt)e'- if it were a close neighbor alignment we would 
expect that (jt)e (the average tangential ellipticity for all 
the pair members in each pair separation bin 9) carries 
all of the signal, which is not the case. In fact we find 
(lt)e ~ 0-2(7t7t)e, which means that close neighbor effect 
can hardly exceed 20% of the small scale signal. 
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Fig. 14. Tangential (top panel) and radial (bottom panel) 
shear correlation functions (7*74)0 and (7 r 7 r }# down to 
3". The solid, long-dashed (hidden by the solid line) and 
short-dashed lines are predictions from the same models 
as in Figure @. 



A forthcoming paper using the same data set will be 
devoted to the measurement of E and B modes (as defined 
in Crittenden ct al. 2000a), and we will study this small 



scale signal in more detail. At this stage of the analysis we 
cannot exclude a possible residual systematic. However, a 
preliminary analysis shows that the B mode down to 3" 
is much smaller than the E mode, which is hard to have 
if the signal comes from residual systematics. 

8. Conclusion 

Using 6.5 sq. deg. of the VIRMOS survey in progress at 
the CFHT, we were able to measure various 2-points cor- 
relation statistics of cosmic shear. The top-hat variance, 
the aperture mass statistic and different shear correla- 
tion functions gave consistent results over a wide range 
of scales. Further tests of the lensing origin of the signal 
exploiting the scalar nature of the gravitational potential 
were also convincingly verified. We demonstrated that the 
level of systematics, in particular the intrinsic alignment 
of galaxies, is likely to be small, and does not contribute to 
the signal for scales larger than 1'. We believe that these 
results demonstrate the significance of our detection of 
shear correlations due to gravitational lensing. The qual- 
ity of the data and the adequate size of our survey allow 
us to constrain cosmological models of the large-scale dis- 
tribution of dark matter in the universe. 

We have obtained tight constraints on the cosmologi- 
cal parameters VLq and erg . These results suggest that high 



12 



Van Waerbeke et al: Cosmic Shear Statistics and Cosmology 




0.005 

5 10 15 20 

0[arcsec] 

Fig. 15. Same as Figure |lj but for the shear correlation 
function (77)9. 



precision measurements can be made with larger surveys 
on a much larger set of cosmological parameters. The final 
stage of the VIRMOS survey is to accomplish 16 sq. deg. 
in patches of 4 sq. deg., 4 colors each, thus allowing the 
possibility to use the photometric redshifts of the galax- 
ies. The use of photometric redshift will not only improve 
the scientific interpretation of cosmic shear (e.g. doing to- 
mography as in |Hu 1999| ) but will be useful to measure the 
intrinsic alignment itself (which can be used to constrain 
galaxy formation models for instance). 

The constraints obtained so far are within a framework 
of structure formation through gravitational instability 
with Gaussian initial conditions and Cold Dark Matter. 
As the amount of observations increases and the measure- 
ment quality improves, the first hints of the shape of the 
power spectrum will be soon available. It opens new means 
of really testing the formation mechanisms of the large- 
scale structure and the cosmological parameters beyond 
the standard model (Uzan fc Bernardeau 2000). 



Over the last two years, we have seen the transition 
from the detection of the weak lensing signal to the first 
measurements of cosmological parameters from it. The 
agreement between theoretical predictions and observa- 
tional results with such a high precision indicates that 
the measurement of cosmic shear statistics is becoming 
a mature cosmological tool. Many surveys are under way 
or scheduled for the next 5 years. They will use larger 
panoramic cameras than the CFH12K, and will cover solid 
angles 10 to 100 times wider than this work. The results of 
this work give us confidence that cosmic shear statistics 
will provide valuable measurements of cosmological pa- 
rameters, probe the biasing of mass/light, produce maps 



of the dark matter distribution and reconstruct its power 
spectrum. 
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